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Abstract 

A measurement of the inclusive bottom jet cross section is presented for events containing a Z 
boson in pp collisions at y/s = 1.96 TeV using the Collider Detector at Fermilab. Z bosons are 
identified in their electron and muon decay modes, and b jets with Et > 20 GeV and \n\ < 1.5 are 
identified by reconstructing a secondary decay vertex. The measurement is based on an integrated 
luminosity of about 330 pb _1 . A cross section times branching ratio of a(Z + b jets) x B(Z — ► 
£ + £~) = 0.93 ± 0.36 pb is found, where B(Z — > £ + £~) is the branching ratio of the Z boson or 
7* into a single flavor dilepton pair (e or fi) in the mass range between 66 and 116 GeV/c 2 . The 
ratio of b jets to the total number of jets of any flavor in the Z sample, within the same kinematic 
range as the b jets, is 2.36 ± 0.92%. Here, the uncertainties are the quadratic sum of statistical 
and systematic uncertainties. Predictions made with NLO QCD agree, within experimental and 
theoretical uncertainties, with these measurements. 
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I. INTRODUCTION 



The measurement of the Z + b jet production cross section provides an important test 
of quantum-chromodynamics (QCD) calculations [1]. The cross section is sensitive to the b 
quark density in the proton and thus tests the perturbative calculations of this quantity. A 
precise knowledge of the b quark density is essential to accurately predict the production of 
particles that couple strongly to b quarks including Higgs bosons (h) within supersymmetry 
models (gb — > hb, bb — > h) [2, 3] or single top production [4] within the standard model 
(qb — > q't and gb — > Wt). The Z + b jet cross section is also an important test of the 
background predictions to standard model Higgs boson production in association with a Z 
boson, ZH -> Zbb [5]. 

The Z cross section [6-8] and Z+jets cross section [9] have been measured at the Tevatron. 
Next to leading order (NLO) QCD calculations are found to describe the data. In this paper, 
the first measurement of the b jet cross section for events with a Z boson using the Collider 
Detector at Fermilab (CDF) [10] is reported. A similar measurement has been made recently 
by the D0 collaboration [11]. 

The dominant production diagrams are gb — > Zb and qq — > Zbb: in NLO calculations 
they contribute about 65% and 35%, respectively. At present the b quark density is derived 
from the gluon distribution function [12] and agrees well with the available measurements 
of the contribution to the proton structure function F 2 for Q 2 < 1, 000 GeV 2 [13], where Q 2 
is the momentum transfer squared. The measurement reported in this paper is sensitive to 
parton densities at higher values with Q 2 approximately equal to the square of the Z mass 
(M z ). 

The analysis uses Run II pp collision data from CDF taken up to September 2004 at a 
center of mass energy of ^/s = 1.96 TeV. The measurement is made by searching for pairs 
of electrons or muons with an invariant mass consistent with Mz and jets which contain a 
displaced secondary vertex consistent with the decay of a long-lived bottom hadron. The 
light (u, d, s, and gluon) and charm (c) jets remaining after this vertex requirement are 
distinguished from the b jets using the mass distribution of the charged particles forming 
the secondary vertex. This technique exploits the larger mass of the b quark compared 
with light and c quarks. The Z cross section is defined to include the irreducible Drell-Yan 
contribution 7* — > within the dilepton invariant mass range 66 < Mu < 116 GeV/c 2 . 
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Note that this cross section is numerically only 0.4% higher than the inclusive Z cross 
section independent of the mass range [6]. The Z + b jet production cross section is defined 
to be proportional to the number of b jets with jet tranverse energy > 20 GeV and 
pseudorapidity \rf ot \ < 1.5 contained in events with a Z boson. 

In Sec. II a brief description of the CDF detector is given, and in Sec. Ill the Monte 
Carlo simulation is described. Section IV summarizes the event selection and the background 
sources. In Sec. V the fraction of b jets within the data sample is determined. In Sec. VI the 
method to measure the cross section is described, and in Sec. VII the sources of systematic 
uncertainties are discussed. The results of the measurement are given in Sec. VIII, and a 
conclusion is presented in Section IX. 

II. THE CDF II DETECTOR 

The CDF II detector is described in detail elsewhere [10]. It is a general purpose, nearly 
hermetic detector situated around the pp collision point. A coordinate system is used, in 
which 9 is the polar angle with respect to the proton beam direction, is the azimuthal 
angle and i] = — lntan(6 l /2) is the pseudorapidity. The transverse energy and transverse 
momentum of a particle is defined as E? = E sin 9 and Pt = P sin 9, respectively, where E is 
the energy measured by the calorimeter and p is the momentum measured in the tracking 
system. The missing transverse energy vector is defined as $ T = — J2i E^rii, where rfj is a 
unit vector that points from the interaction vertex to the center of the ith calorimeter tower 
in the transverse plane and E l T is the transverse energy of the ith tower. The quantity $ T 

— * 

is the magnitude of $ T , which is corrected for all identified muons in an event [7]. 

The transverse momenta of charged particles are measured by an eight-layer silicon strip 
detector [14-16] and the central outer tracker (COT), a 96-layer drift chamber [17] located 
inside a solenoid that provides a 1.4 T magnetic field. The innermost layer of the silicon 
detector is located on the beryllium beampipe at a radius of 1.5 cm, and the outermost 
layer is located at 28 cm. The silicon detector provides tracking in the pseudorapidity 
region |^| < 2, with partial coverage up to \r}\ < 2.8. The single hit resolution is about 
11 fim. Located outside of the silicon detector, the COT is a 3.1 m long, open-cell drift 
chamber with an active tracking region extending radially from 41 cm to 137 cm. The COT 
provides coverage for \r)\ < 1. For tracks with px > 1.5 GeV/c and silicon hits the resolution 
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on the impact parameter is about 34 jum [18], including the transverse size of the beam of 
about 25 (Urn. 

Located outside the solenoid, a segmented sampling calorimeter is used for the mea- 
surement of particle energies. The central part of the calorimeter covers the region 
< 1.1 [19, 20], and the forward part of the calorimeter consists of two identical de- 
tectors covering 1.1 < \rj\ < 3.6 [21]. The central calorimeter uses lead-scintillator sampling 
in the electromagnetic compartment and steel-scintillator sampling in the hadronic compart- 
ment. It is instrumented with proportional strip and wire chambers [central electromagnetic 
shower maximum detector (CES)]. They are located at a depth of about six radiation lengths 
where the lateral profile of electromagnetic showers is expected to be maximal and have a 
segmentation of 1.5 cm. The forward calorimeter uses lead-scintillator sampling for the 
electromagnetic compartment and iron-scintillator for the hadronic compartment. Further 
details about the calorimeters can be found in Ref . [22] . 

Drift chambers, located outside the central hadron calorimeters and behind 60 cm of 
iron shielding, detect muons with \r]\ < 0.6 [23]. Additional drift chambers and scintillation 
counters detect muons in the regions 0.6 < \r]\ < 1.0 and 1.0 < \r]\ < 1.5. Gas Cherenkov 
counters [24] measure the average number of pp inelastic collisions per bunch crossing to 
determine the luminosity. 

III. MONTE CARLO SIMULATION 

A Monte Carlo simulation is used to correct for inefficiencies due to the selection re- 
quirements and detector effects. The Monte Carlo generator pythia v6.2 [25] is used to 
generate the Drell-Yan signal and the background processes, using CTEQ5L parton density 
functions (PDFs) [26]. An underlying event model that describes the interactions of the 
spectator partons and initial state QCD radiation has been included in the generation. This 
model has been tuned to describe the Tevatron data [22, 27]. The decays of the b hadrons 
are generated by the Monte Carlo generator QQ v9.1 [28]. The CDF detector response is 
simulated using a geant based detailed detector simulation [29, 30]. 

The Drell-Yan Monte Carlo samples are normalized to the next to next to leading order 
(NNLO) QCD cross section of 251.3 pb [31] for 66 < M u < 116 GeV/c 2 , so that comparisons 
with the data can be made. For ti and ZZ processes, which contribute to the background, 
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the NLO QCD cross sections are used for the normalization: a t t = 6.77 pb [32] and a Z z = 
1.4 pb [33]. 

Simulated events are reconstructed in the same manner as the data events, and the same 
event selection criteria are applied. 

IV. EVENT SELECTION AND BACKGROUND 
A. Event Selection 

Z bosons are detected in their decays into two electrons or two muons with an invari- 
ant mass of the two leptons M u (I = e or /i) between 66 and 116 GeV/c 2 . The trigger 
requirements and the lepton selection follow closely those described in detail in Ref. [7]. 

Electrons are triggered by requiring a cluster of electromagnetic energy with Et > 18 GeV 
and 1 77 1 < 1.1 matched to a track with px > 10 GeV/c. Further requirements are made on 
position matching and the shower shape in the CES. At least one trigger electron candidate is 
required. The second electron candidate can either be in the central or forward calorimeter, 
and looser identification criteria are imposed. For forward electron candidates no matching 
track is required. All electron candidates are required to be isolated from other calorimeter 
energy deposits [7]. 

Muons are triggered by requiring a track with px > 18 GeV/c and < 1.0 and a 
track segment in the muon chambers that matches the extrapolated position of the track. 
At least one muon candidate that satisfies the trigger requirements is required. The other 
muon candidate is not required to have signals in the muon chambers. All muon candidates 
are required to have a calorimeter energy deposit consistent with that of a muon and to 
be isolated from other energy depositions [6, 7]. The two highest px muon candidates are 
required to have opposite electric charges. 

Candidate Z boson events are selected if 66 < Mu < 116 GeV/c 2 . A total of 27,659 
candidate events are observed in the electron channel and 15,698 events in the muon channel. 

Having selected an event with a Z boson candidate, jets with E^ > 20 GeV and 
Irf^l < 1.5 are searched for. Jets are defined by a cone jet algorithm with a cone size 
R = a/ (At/) 2 + (A0) 2 = 0.7 [22]. The jet energy is corrected to the hadron level energy. 
The hadron level energy is defined to include all particles from the pp collision within the 
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jet cone, including particles from the hard scatter, multiple parton-parton interactions, and 
beam remnants. The jet energy is also corrected for particles produced in additional pp 
interactions, reconstructed in the same bunch crossing. The jets are not corrected to the 
parton level to be independent of the Monte Carlo modeling of this correction. 

Events must satisfy either ~$ T < 25 GeV or H T < 150 GeV, where H T is the scalar sum 
of $ T and the transverse energies of all leptons and jets in the event [34]. This requirement 
reduces background from tt events in which the t decays to Wb and both the W bosons decay 
leptonically (W — > Iv) by about 80%, while reducing the signal by only 4%, as determined 
from Monte Carlo simulation. 

A b jet is defined as any jet that has at least one b hadron within a cone of 0.7 around 
the jet axis. In this analysis a b jet is identified through the presence of a displaced vertex 
within the jet arising from the decay of the long-lived bottom hadron. The algorithm used 
was optimized for the measurement of the top quark production cross section [35] but found 
to give adequate efficiency and purity for the present analysis. A jet that has a reconstructed 
displaced vertex is called a "6 tagged" jet. The displaced vertex algorithm uses a two-pass 
approach to find a secondary vertex. In the first pass an attempt is made to reconstruct a 
secondary vertex using one track with px > 1-0 GeV/c and two or more additional tracks 
with px > 0.5 GeV/c, and all tracks are required to have an impact parameter significance 
do/o'do > 2-5. Here, d is the minimum distance between the track and the primary vertex 
in the plane transverse to the beam direction and has uncertainty <Jd - If the first pass is 
unsuccessful, a second pass is made using two tracks with px > 1.5 GeV/c for one track, 
Pt > 1.0 GeV/c for the other and d /a d() > 3 for both. The jet is labeled as a tag if 
the transverse displacement significance \L 2 d\/o'l 2D > 7.5. Here, L 2 d is the distance from 
the primary vertex to the secondary vertex in the plane transverse to the beam direction 
projected onto the jet axis, and <7l 2D i s the estimated uncertainty. The distance L 2D is 
defined as positive if the angle between the transverse displacement and the jet direction is 
less than 90°, and as negative otherwise. A positive tag has positive L 2 d, and a negative tag 
has negative L 2 d- 

The b tagging efficiency is determined as a function of and r]i ct from a separate data 
sample of about 28,000 dijet events, where one of the jets has a reconstructed semileptonic 
b or c decay [35]. The ratio of the efficiency in data to that in Monte Carlo simulation is 
found to be 0.91 ± 0.06. The average data b tagging efficiency for jets in this analysis is 
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33 ± 2%. 

In total, 115 tagged jets are selected, 60 in the electron channel and 55 in the muon 
channel. This compares with a Monte Carlo estimate of 69 in the electron and 45 in the 
muon channel. The Monte Carlo estimate includes the Drell-Yan contribution, which has 
been scaled by the factors obtained in the fit to the secondary vertex mass (see Section V), 
and the backgrounds listed in Section IV B. Out of the 115 tagged jets 16 are negatively 
tagged compared with a Monte Carlo estimate of 16.7. One event contains two positively 
tagged b jets compared with a Monte Carlo expectation of 1.46. 

B. Backgrounds 

Backgrounds to Z + b production can arise from misidentified leptons, from genuine lep- 
tons and b jets coming from other processes, or from light jets or c jets that are misidentified 
as b jets. The first two background sources are discussed in this section, and the latter is 
discussed in section V. 

The background in which one or both reconstructed electrons in the Z — > e + e~ channel 
are misidentified from other particles in the final state is estimated from the data. The 
probability that a jet will pass all electron identification criteria is determined from several 
inclusive jet samples. These samples have negligible prompt electron content. The proba- 
bilities are parameterized as a function of and are on average 0.1% for central jets and 
1.5% for forward jets. Because an associated track is not required for a forward electron, 
the misidentification probability is much higher than for central electrons. 

A sample of events in which exactly one trigger electron is reconstructed is now taken. 
The trigger electron is paired with any other jet in the event such that the invariant mass 
of the electron and the jet lies within the Z mass window. The jet energy is taken at 
the electromagnetic calorimeter energy scale i.e. the correct energy scale for an electron or 
photon. A weight, which equals the jet misidentification probability, is assigned to the jet. 
The total background to inclusive Z production is then the sum of all weights. If there is 
more than one jet in the event that forms an invariant mass within the mass window of the 
Z, each combination is used. Background distributions are derived by weighting the electron 
+ jets distributions with the weights. The background for Z+jets (Z + b jets) is the sum of 
those weights for events which contain at least one jet {b jet) in addition to the one paired 
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with the electron. 

Using this method, a background contribution of 3.1 ± 1.5 b tagged jets is estimated 
within the Z mass window. 

Background in the Z — > channel in which the electric charges of the two muon can- 

didates are uncorrelated is estimated from events with two reconstructed muon candidates 
that have the same electric charge [7] . This background comes from events in which one or 
both muon candidates arise from hadron decays or events in which hadrons in the final state 
are misidentified as muons. Due to the low statistics of the events with a b tagged jet, the 
fractional contribution of this background is estimated as the observed ratio of generic jets 
(i.e., any jet regardless of which quark flavor or gluon it originated from) with a like-sign Z 
candidate to those with an unlike-sign Z candidate. The number of b tagged jets from this 
background source is estimated to be 0.24 ± 0.12. This number is in agreement with the 
zero tagged jets with a like sign Z candidate observed in the data. 

The background from other processes is estimated from Monte Carlo simulations. The 
production of ti pairs is found to contribute 0.25±0.05 (0.24±0.05) tagged jets to the electron 
(muon) channel. The production of ZZ is found to contribute 0.36±0.07 (0.28±0.06) tagged 
jets to the electron (muon) channel. Backgrounds from other processes such as WW, WZ, 
bb production or Z — > t + t - are estimated to be negligible. 

The invariant mass of the dilepton pair is shown in Figure 1 for events with at least 
one generic jet with > 20 GeV and |r] jct | < 1.5. The data are compared with the 
expected Drell-Yan and estimated background contributions. Good agreement is observed 
near M z and in the tails. Figure 2 shows the dilepton invariant mass for events with at least 
one positively tagged jet in the same kinematic range. The data are well modeled by the 
simulation, and a clear Z signal is also observed in this sample. 

V. FRACTION OF b JETS 

The fraction of b jets in the tagged jet sample is estimated by performing a fit to the 
invariant mass of all charged tracks attached to the secondary vertex M s [36]. On average 
b jets have a larger M$ than c jets or light jets due to the larger mass of b hadrons. The 
Ms distribution for both positively and negatively tagged jets is used, in order to better 
discriminate between heavy and light jets, since jets with a genuine secondary vertex over- 
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Figure 1: The invariant mass of the dilepton pair for the sample with jets with E£ > 20 GeV 
and \rf et \ < 1.5 compared with the expectation from signal and background sources. The Drell- 
Yan Monte Carlo has been normalized to the luminosity of the data sample assuming the NNLO 
Drell-Yan cross section. 

whelmingly have positive tags, whereas jets with a false secondary vertex may have positive 
or negative tags. In addition using the negative tags in the fit allows a better separation 
of the charm and light quark contributions that have similar Ms distribution for positively 
tagged jets. 

The Ms distributions for positively and negatively tagged data jets are shown in Figure 3. 
The distributions of Ms for b, c, and light jet events are taken from the Monte Carlo 
simulation. The Monte Carlo light jet distribution has been corrected using data as described 
below. 

It is important to ensure that the simulation models the Ms distributions well. The Ms 
distributions are affected by the tracking efficiency, the charged particle multiplicity of the 
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Figure 2: The invariant mass of the dilepton pair for the sample with positively tagged jets with 
Ej?" > 20 GeV and \rf et \ < 1.5 compared with the expectation from signal and background sources. 
The Drell-Yan Monte Carlo has been scaled by the factors determined in Section V. 



b hadron decay, and the fraction of jets with two b or c hadrons. Systematic uncertainties 
have been assigned to each of these contributions (see Section VII), and the Ms distribution 
in a Monte Carlo simulation of b jets has been compared with the dijet sample used to 
determined the b efficiency (see Section IV). The simulation is found to describe the data 
within the uncertainties quoted [37]. 

The number of positive tags in light quark jets is larger than the number of negative 
tags due to long lived particles, such as Kg and A, and due to vertices produced by nuclear 
interactions of particles with material in their path [35]. The current simulation may not 
describe these effects accurately. In order to investigate these effects, a sample of approx- 
imately 13,000 7 + jet data events is taken, with jets in the same kinematic range as for 
the Z + b jets selection. This sample was chosen since it has higher statistics than and a 
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Figure 3: The mass at the secondary vertex, Mg, for (a) positively and (b) negatively tagged jets 
with > 20 GeV and \rf et \ < 1.5. The non-Drell-Yan background has been subtracted from 
the data. The data are compared with the sum of the light, c and b Monte Carlo templates after 
being scaled by the factors pi, p c , and pb, respectively. The open white area represents the light 
quark template, the lightly shaded are the c quark template and the dark shaded are the b quark 
template. 
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similar event topology as the Z + b jets sample. The M s distribution of the 7 + jet data 
is obtained separately for positive and negative tags. The fraction of light, c, and b jets in 
each sample is determined using Monte Carlo distributions as templates and performing a 
likelihood fit, similar to the one described below, in which the normalization of each Monte 
Carlo template is allowed to vary. It is found that the ratio of positive to negative tags 
for light jets in the data is 1.49, compared with the expectation from the Monte Carlo of 
1.93. To correct for the difference between the data and the Monte Carlo simulation, the 
Monte Carlo light negative tagged jets are increased by a factor of 1.93/1.49 = 1.30. This 
reweighted template shape is taken for the central value, and the full difference from unity, 
±30%, is taken as the systematic uncertainty. 

A binned maximum likelihood fit using Poisson statistics is performed to the Ms dis- 
tributions of the positively and negatively tagged jets in the Z sample for range Ms < 
3.5(1) GeV/c 2 for positively (negatively) tagged jets. The range is chosen in order to have 
enough statistics for the fit. The Monte Carlo distributions are taken as templates and scaled 
by factors pb, p c , and pi for b, c, and light jets, respectively. The data distribution is used 
after subtraction of the non-Drell-Yan background as estimated in Section IV B. The quan- 
tity pi, is thus the number of fitted reconstructed signal b jets in data divided by the number 
in the simulation. Only statistical errors from the data and Monte Carlo are used in the fit. 
The fit takes into account the Monte Carlo statistical errors using the method described in 
Ref. [38]. The fit gives values of p b = 0.93 ± 0.29, p c = 1.69 ± 0.94 and p x = 1.36 ± 0.53. The 
correlation coefficient between pb and p c is —0.68 and between pb and pi is 0.10. 

The number of b jets in the Z sample after subtracting the background contributions is 
estimated to be N-£, ata (Z + b jets) = 45 ± 14. A check is performed by fitting only positively 
tagged jets, and good agreement is obtained with a value of pb = 0.95 ± 0.31. 

VI. CROSS SECTION 

The inclusive b jet cross section a{Z + b jets) is proportional to the number of b jets with 
> 20 GeV and |?7 jet | < 1.5 and is defined for dilepton masses 66 < M u < 116 GeV/c 2 . 
The branching fraction B(Z — > is defined for a single lepton flavor. All the cross 

sections and cross section ratios presented are fully corrected for detector response and 
resolution, and presented at the hadron level. No corrections are made for the underlying 
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event or for hadronization effects (see Section IV). They are defined for jets with a cone size 
of 0.7 that include a b hadron within this cone. 

A ratio method is used to extract the cross section. In doing so use can be made of the 
uncertainties estimated for the inclusive Z cross section measurement [7] for the lepton and 
trigger selection. First a measurement is made of the ratio of the Z + b jet cross section to 
the total Z cross section: 

<x(Z + bjets) _ N Bata (Z + 6jets)/e(Z + 6 jets) 

<r(Z) N Data (Z)/e(Z) ' 1 j 

where Nn ata (Z + 6 jets) is the fitted number of b jets (see Section V) and Nv ata (Z) is the 

total number of events with a lepton pair in the mass range 66 < Mm < 116 GeV/c 2 in the 

data. In both cases the number of data events is taken after subtraction of the background 

contributions (see Section IV B). The efficiencies of the Z + b jet and the Z boson selections 

are e(Z + 6 jet) = 7.7% and e(Z) = 27%, respectively. They are determined from pythia 

Monte Carlo simulation and are corrected for any differences from the data. The ratio 

e(Z + b jets) / e(Z) is also determined using HERWIG v6.5 [39] and a similar result is obtained 

(0.286 for herwig compared with 0.285 for pythia). 

The cross section is then calculated as: 

<t(Z + fcjets) x B(Z -> ttr) = g ( Z +y s ) . acDF ( Z ) x B (Z -> £+r), (2) 

where a CJiF (Z) x B{Z -> £+£') = 254.9 ± 3.3(stat.) ± 4.6(syst.) ± 15.2(lum.) pb is the CDF 
measurement of the Z production cross section times branching fraction for a single lepton 
flavor [6]. 

A measurement is also made of the ratio of the Z + b jets cross section to the Z+ 
generic jets cross section in order to measure the fraction of jets that contain at least one b 
hadron. The Z+ generic jets cross section is proportional to the number of generic jets with 
> 20 GeV and |r] jct | < 1.5. The ratio of the Z + b jets to Z+ generic jets cross section 
is obtained from the data as: 



g(Z + b jets) = N Data (Z + b jets)/ e(Z + b jets) 

<7(Z+jets) iV Data (Z+jets)/e(Z + jets) ' U 
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where Nn ata (Z + jets) is the number of generic data jets for events with a lepton pair within 
the Z mass window after background subtraction and e(Z + jets) = 24% is the efficiency of 
the Z + jet selection. 

The validity of the kinematic distributions in the Monte Carlo simulation was checked by 
comparing the transverse energy and pseudorapidity distributions for generic jets and for 
the positively b tagged jets. These distributions, which are shown in Figure 4 for generic 
jets and in Figure 5 for b jets, demonstrate that pythia describes the data well and thus 
may be used to correct the data. 

VII. SYSTEMATIC UNCERTAINTIES 

The total systematic uncertainty on the measurement is estimated by adding the uncer- 
tainties described in this section in quadrature. 

The uncertainty on the jet energy scale was estimated by changing the scale in the Monte 
Carlo simulation by the uncertainty following the procedure of Ref. [22]. An additional 
uncertainty of ±0.6% [40] is applied on the b jets to account for possible differences in the 
description of b jets by the Monte Carlo simulation compared with light jets. 

Possible differences between the Monte Carlo simulation and data E 3 ^ 1 and rf et distribu- 
tions, which may arise for many reasons including uncertainties in the parton distribution 
functions, inadequacies in the QCD model or uncertainties in the b quark fragmentation, 
are accounted for in the systematic error by weighting the Monte Carlo hadron level distri- 
butions by (E 3 ^l Lad /35) ±aE T and (IVhLiI + 0-5) ±a,? , respectively. Here -E^ Had is the transverse 
energy in GeV, and is the pseudorapidity of the Monte Carlo hadronic jet. The values 
of 35 GeV and 0.5 units of pseudorapidity are arbitrary and chosen to be close to the center 
of the distributions. The parameter otE T is set to 0.2 (1.0) for untagged (tagged) jets. The 
parameter a v is set to 0.2 (1.0) for untagged (tagged) jets. These parameters are chosen 
by comparing the weighted and unweighted Monte Carlo reconstructed distributions with 
those of the data and choosing the weight such that it is still statistically compatible with 
the data. 

The b jet tagging efficiency has a 6.6% uncertainty [35]. 

Since any changes in the shapes of the light, c and b jet templates can change the value of 
Pb obtained, various possible variations of the template shapes are included in the systematic 
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Figure 4: The and rf et distributions for generic jets with E8? > 20 GeV and \rf et \ < 1.5. The 
Drell-Yan Monte Carlo has been scaled such that the total number of jets in the simulation is the 
same as in the data. 2^ 
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uncertainty. An uncertainty in the negative to positive light jet tag rate is estimated by 
weighting the negative tags of the light template by ±30%. This is the observed difference 
between data and Monte Carlo simulation for the negative to positive light jet tag ratio 
(see section V). Since it is not a priori known whether the jets contain one or two heavy 
quarks, an uncertainty in the shape of the c (6) quark template is estimated by taking 
half the difference between results obtained using the templates when there is one or two c 
(b) hadrons within a cone of 0.7 around the jet axis. An uncertainty due to possible track 
reconstruction inefficiencies is estimated by recalculating templates after randomly removing 
3% of the tracks in the Monte Carlo simulation. An uncertainty in the shape of the b 
quark template due to the uncertainty in the mean b hadron charged particle multiplicity is 
estimated by weighting the Monte Carlo multiplicity such that the mean changes by ±0.15, 
which is the difference between the Monte Carlo simulation used in this analysis and data 
measurements [41]. This systematic uncertainty also covers any differences observed between 
data and simulation in the 6-quark rich sample used to determine the efficiency (see Section 
IV). 

An uncertainty on the misidentified electron and muon backgrounds is estimated to be 
±50% of the estimated contribution (see Section IV B). An uncertainty in the ZZ and ti 
backgrounds of ±20% is estimated. These uncertainties include the NLO cross section uncer- 
tainty estimated by varying the factorization scale by a factor two [32, 33], the cross section 
uncertainty due to the uncertainty on the top quark mass, and experimental uncertainties 
(e.g. the lepton identification and b jet tagging efficiencies). 

A 1.8% uncertainty on the selection efficiency of Z events is taken from Ref. [6]. A 5.8% 
uncertainty is taken for the luminosity determination [42]. 

The effect of these contributions on the systematic uncertainty on the cross section is 
listed in Table I. The systematic uncertainties are evaluated separately for the cross section 
ratio measurements and are also included in Table I. The total systematic uncertainty on the 
cross section is 22% compared with a statistical precision of 31%. The largest uncertainty 
arises from the assumed Monte Carlo distribution and will be reduced with higher data 
statistics to constrain this distribution. 
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Table I: The systematic uncertainties on the cross section and ratio measurements. The total 
systematic uncertainty on each measurement is estimated by adding the individual uncertainties 
in quadrature. 



Source of Uncertainty 


S(a(Z + 6 jets)) 
(%) 


5(a(Z + b jets) /a (Z)) 
(%) 


5(a(Z + b jets) /a(Z + jets)) 
(%) 


jet energy scale 


3.9 


3.9 


6.7 


b jet energy scale 


0.2 


0.2 


0.2 


MC rf et dependence 


6.6 


6.6 


6.6 


MC dependence 


16.1 


16.1 


16.2 


b tagging efficiency 


6.6 


6.6 


6.6 


light jet template 


1.8 


1.8 


1.8 


single/double c quark in jet 


2.8 


2.8 


2.8 


single/double b quark in jet 


3.6 


3.6 


3.6 


track reconstruction efficiency 


9.9 


9.9 


9.9 


b hadron multiplicity 


1.0 


1.0 


1.0 


misidentified lepton background 


1.3 


0.9 


0.4 


other backgrounds 


0.3 


0.3 


0.3 


Z selection efficiency 


1.8 


n. a. 


n. a. 


luminosity 


5.8 


n. a. 


n. a. 


total 


22.9 


22.1 


22.8 



VIII. RESULTS 

The b jet cross section for events with 66 < M« < 116 GeV/c 2 and E 1 ^ > 20 GeV and 
\rf et \ < 1.5 is measured to be 

<j(Z + b jets) x B{Z -> tT) = 0.93 ± 0.29(stat.) ± 0.21(syst.) pb. 

The theoretical prediction of 0.45 ± 0.07 pb is consistent with the measurement. The theo- 
retical cross section and the ratios listed below were evaluated at NLO using MCFM [1, 43] 
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with the parton distribution functions from CTEQ6M [44]. The factorization and renor- 
malization scales were set to M z . The uncertainty on the NLO prediction includes the 
uncertainty arising from the choice of these scales, which were varied from M z /2 to 2M Z: 
the uncertainty on the parton distribution functions and the uncertainty on the strong cou- 
pling constant. 

The ratio of Z + fejets production to inclusive Z production is measured to be 

a{Z + 6jets)/<r(Z) = 0.0037 ± 0.0011(stat.) ± 0.0008(syst.). 

The pythia estimate of 0.0035 and the NLO QCD calculation of 0.0019 ±0.0003 both agree 
with the measured value. 

The NLO QCD cross section a(Z + b jets) and the ratio a(Z + b jets) /a(Z) do not include 
the effects of underlying event and hadronization. These effects are estimated from pythia 
to change the NLO QCD cross section and the ratio to the inclusive Z cross section by 
+ 10% for the underlying event and —1% for the hadronization. 

The ratio of Z + frjets production compared to Z + jets production is measured to be 

a(Z + b jets) /a(Z + jets) = 0.0236 ± 0.0074(stat.) ± 0.0053(syst.) 

The estimate from pythia of 0.0218 and the NLO QCD calculation of 0.0181 ± 0.0027 
agree with this result. The NLO QCD cross section has not been corrected for underlying 
event and hadronization. These effects are estimated from pythia and change the value 
by —7% and +1%, respectively. The DO collaboration [11] obtained a measurement of 
0.021 ± 0.004(stat.)^QQQ3(syst.) for the fraction of events with at least one b jet, with P^ et > 
20 GeV/c and \rf ct \ < 2.5, to those events with at least one generic jet within the same 
kinematic range. Note that in the analysis by the DO collaboration, an NLO QCD prediction 
for the ratio of the Z+b quark to the Z+c quark production rate is assumed. This assumption 
results in smaller uncertainties for the DO analysis despite the present analysis having more 
than twice the number of tagged b jets and similar b purity. The present analysis does not 
take the approach followed by the DO collaboration in order to maintain independence of 
the measurement from QCD calculations. 
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IX. CONCLUSIONS 



In conclusion, the inclusive 6 jets production cross section in events with a Z boson has 
been measured for jets with E T > 20 GeV and \rj^ et \ < 1.5 in pp collisions at y/s = 1.96 TeV 
at the Tevatron using CDF data with an integrated luminosity of about 330 pb _1 . 

The Z + b jets cross section times branching ratio is measured as 0.93 ± 0.36 pb for the Z 
boson decaying into a single charged lepton flavor. The ratio of the Z + frjets to the Z + jets 
cross section is 2.36 ± 0.92%. This is the first measurement of the Z + frjets cross section 
without any assumptions on the Z + cjets production rate. The uncertainty is dominated 
by the limited statistical precision, which is expected to improve significantly in the near 
future. 

Within current uncertainties, the theoretical predictions agree with these measurements, 
which are sensitive to the b quark density at high momentum transfer. 
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